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ABSTRACT 

We extend the 3-point intrinsic alignment self-calibration technique to the gravita- 
tional shear-intrinsic ellipticity-intrinsic ellipticity (Gil) bispcctrum. The proposed 
technique will allow the measurement and removal of the Gil intrinsic alignment 
contamination from the cross-correlation weak lcnsing signal. While significantly de- 
creased from using cross-correlations instead of auto-correlation in a single photo-z 
bin, the Gil contamination persists in adjacent photo-z bins and must be accounted 
for and removed from the lensing signal. We relate the Gil and galaxy density-intrinsic 
ellipticity-intrinsic ellipticity (gll) bispectra through use of the galaxy bias, and de- 
velop the estimator necessary to isolate the gll bispectrum from observations. We 
find that the Gil self-calibration technique performs at a level comparable to that of 
the gravitational shear-gravitational shear-intrinsic ellipticity correlation (GGI) self- 
calibration technique, with measurement error introduced through the gll estimator 
generally negligible when compared to minimum survey error. The accuracy of the re- 
lationship between the Gil and gll bispectra typically allows the Gil self-calibration 
to reduce the Gil contamination by a factor of 10 or more for all adjacent photo-z bin 
combinations at £ > 300. For larger scales, we find that the Gil contamination can be 
reduced by a factor of 3-5 or more. The Gil self-calibration technique is complemen- 
tary to the existing GGI self-calibration technique, which together will allow the total 
intrinsic alignment cross-correlation signal in 3-point weak lensing to be measured and 
removed. 
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1 INTRODUCTION 

Weak gravitational lensing due to large scale structure (cosmic shear) has become a promising source of cosmological in- 
formation. A new generation of ground- and space-based surveys suited for precision weak lensing measurements have been 
developed with the importance of this new probe in mind. These ongoing, future, and proposed surveys (e.g. CFHTLS 1 , DES 2 , 
EUCLID 3 , HSC 4 , HST 5 , JWST 6 , LSST 7 , Pan-STARRS 8 , and WFIRST 9 ) promise to provide greatly improved measurements 
of cosmic shear using the shapes of up to billions of galaxies. There has been much work done to explore the potential of these 
cosmic shear measurements, which we review in Troxel & Ishak (2012a), for both the 2- and 3-point cosmic shear correlations. 
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Figure 1. The 3-point galaxy intrinsic alignment correlations. Blue galaxies are intrinsically aligned (I), while red galaxies are lenscd (G). 
The view of the system on the sky is represented in the lower right panels. Each preceding panel demonstrates the galaxy configuration 
at some distinct redshift such that Zi < Zj < z^. 



Beyond the constraints obtained on cosmological parameters from the 2-point cosmic shear correlation and the corre- 
sponding shear power spectrum, the 3-point cosmic shear correlation and shear bispectrum are able to break degeneracies 
between the cosmological parameters that the power spectrum alone does not (Takada & Jain 2003; Vafaei et al. 2010). The 
results of Takada & Jain (2004), for example, showed that the constraints on the dark energy parameters and the matter 
fluctuation amplitude should be able to be improved by a further factor of 2-3 using the bispectrum measured in a deep 
lensing survey. Most recently, parameter constraints were derived by Semboloni et al. (2010) using weak lensing data from 
the HST COSMOS survey, measuring the third order moment of the aperture mass measure. Their independent results were 
consistent with WMAP7 best-fit cosmology and provided an improved constraint when combined with the 2-point correlation. 
The bispectrum also allows us to explore information about non-Gaussianity in the universe that is inaccessible at the 2-point 
level, providing constraints on the degree of non-Gaussianity in addition to the information on other cosmological parameters 
(see for example Matarrese, Verde & Jimenez (2000) ; Verde et al. (2001); Takada & Jain (2004); Jeong & Komatsu (2009); 
Huterer, Komatsu & Shandera (2010); Munshi et al. (2011) and references therein.) 

However, several systematic effects limit the precision of cosmic shear measurements, which must be accounted for in 
order to make full use of the potential of future weak lensing surveys (see for a summary Troxel & Ishak (2012a) and references 
therein), and one of the serious systematic effects of weak lensing is this correlated intrinsic alignment of galaxy ellipticities, 
which act as a nuisance factor (see for example Croft & Metzler (2000); Catelan, Kamionkowski & Blandford (2001); 
Crittenden et al. (2001); Brown et al. (2002); Jing (2002); King & Schneider (2002); Heymans & Heavens (2003); Hirata 
& Seljak (2003b); King & Schneider (2003); Hirata & Seljak (2004); King (2005); Heymans et al. (2006); Mandelbaum 
et al. (2006); Bridle & King (2007); Hirata et al. (2007); Semboloni et al. (2008); Faltenbacher et al. (2009); Okumura 
& Jing (2009); Joachimi & Bridle (2010); Joachimi et al. (2010); Kirk, Bridle & Schneider (2010); Blazek, McQuinn 
& Seljak (2011); Krause & Hirata (2011); Troxel & Ishak (2012b) and references therein). The dark energy equation of 
state can be biased by as much as 50%, for example, if intrinsic alignment is ignored (Bridle & King 2007; Joachimi & 
Bridle 2010). Hirata et al. (2007) found that the matter power spectrum amplitude can be affected by up to 30% due to the 
intrinsic alignment, demonstrating the importance of developing methods to isolate and remove the intrinsic alignment from 
the cosmic shear signal. 

There are two 2-point intrinsic alignment correlations. The first is a correlation between the intrinsic ellipticity of two 
galaxies, known as the intrinsic ellipticity-intrinsic ellipticity (II) correlation, between two intrinsically aligned galaxies. The 
second, known as the gravitational shear-intrinsic ellipticity (GI) correlation, was identified by Hirata & Seljak (2004) and is 
due to a matter structure both causing the alignment of a nearby galaxy and contributing to the lensing signal of a background 
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galaxy. The large-scale GI correlation was first reported by Mandelbaum et al. (2006) in the SDSS, with Hirata et al. (2007) 
finding an even stronger GI correlation for Luminous Red Galaxies (LRGs). The authors showed that this contamination can 
affect the lensing measurement and cosmology by up to 10%, while affecting the matter fluctuation amplitude by up to 30%. 
Heymans et al. (2006) confirmed this through numerical simulations, where a level of contamination of 10% was found. The 
GI correlation was also measured in the SDSS dataset by Faltenbacher et al. (2009) and Okumura & Jing (2009), and most 
recently, Joachimi et al. (2010) measured strong 2-point intrinsic alignment correlations in various SDSS and MegaZ-LRG 
samples. 

Similarly, when we consider the 3-point correlation, the gravitational shear-gravitational shear-gravitational shear (GGG) 
bispectrum suffers from contamination by the 3-point intrinsic alignment correlations. The first correlation is between the 
intrinsic ellipticities of three spatially close galaxies which are intrinsically aligned by a nearby matter structure, known as 
the intrinsic ellipticity-intrinsic ellipticity-intrinsic ellipticity (III) correlation. The second is the gravitational shear-intrinsic 
ellipticity-intrinsic ellipticity (Gil) correlation, where a matter structure which contributes to the lensing of a galaxy in the 
background also intrinsically aligns two spatially nearby galaxies. Finally, there is the gravitational shear-gravitational shear- 
intrinsic ellipticity (GGI) correlation, where a matter structure which intrinsically aligns a galaxy in the foreground also lenses 
two galaxies in the background. The sign of the GGI and Gil correlations can depend both on triangle shape and scale. The 
3-point intrinsic alignment correlations are represented diagrammatically in Fig. I. It was shown by Semboloni et al. (2008) 
that compared to the lensing spectrum, the lensing bispectrum is typically more strongly contaminated by galaxy intrinsic 
alignment, up to 15 — 20% compared to the GGG lensing signal. 

Troxel & Ishak (2012c) demonstrated that the III and Gil intrinsic alignment correlations can be greatly reduced with 
photo-z's by using cross-spectra of galaxies in two different redshift bins so that the galaxies are separated by large enough 
distances to assure that the tidal effect is weak. However, as shown in Fig. 4 of Troxel & Ishak (2012c), while the III bispectrum 
is supressed almost totally for adjacent, large redshift bins, the Gil bispectrum is not. Thus we propose an extension to the 
GGI self-calibration technique (Troxel & Ishak 2012a) for the Gil bispectrum below. There has been significant work in recent 
years toward measuring and removing the effect of galaxy intrinsic alignment on cosmological measurements, and we review 
these briefly. The 2-point intrinsic alignment correlations and methods for their removal have been well studied in recent 
literature (see for example Hirata & Seljak (2004); King & Schneider (2003); King (2005); Bridle & King (2007); Joachimi 
& Bridle (2010); Kirk, Bridle & Schneider (2010); Joachimi & Schneider (2008, 2009, 2010); Zhang (2010a,b)). Using a 
geometrical approach, Shi, Joachimi & Schneider (2010) generalized the nulling technique to the 3-point GGI correlation 
by exploiting its redshift dependence, but it was found that the technique throws out some of the valuable lensing signal. 
Finally, we have recently generalized the approaches of Zhang (2010a) and Zhang (2010b) to develop a set of 3-point self- 
calibration techniques. We proposed using additional galaxy density (cross-)correlations which are already present in weak 
lensing survey measurements to self-calibrate the GGI cross-correlation in Troxel & Ishak (2012a), while the unique redshift 
dependencies of the intrinsic alignment bispectra allowed us to self-calibrate the intrinsic alignment auto-correlations in Troxel 
& Ishak (2012c). 

We propose here an extension to the GGI self-calibration technique of Troxel & Ishak (2012a), developing a technique to 
self-calibrate the Gil cross-correlation bispectrum. While the impact of the Gil cross-correlation on cosmological information 
is significantly less than the impact of the GGI cross-correlation, especially for widely separated redshift bin combinations, 
it still has a contribution for adjacent redshift bin combinations, which must be measured and removed if one hopes to 
achieve precision, systematic free bispectrum measurements. To do this, we develop a relationship between the Gil and galaxy 
density-intrinsic ellipticity-intrinsic ellipticity (gll) bispectrum. We also devise an estimator to measure the gll bispectrum 
from the observed galaxy density-ellipticity-ellipticity bispectrum. We then discuss the performance of such a self-calibration 
technique, describing both measurement and systematic errors due to the self-calibration. 

We organize the paper as follows. In Sec. 2, we briefly discuss the necessary survey parameters and lensing calculations. 
In Sec. 3, we develop the 3-point Gil self-calibration. We first develop a relation between the Gil and gll bispectra, which 
allows for the estimation and removal of the Gil intrinsic alignment cross-correlation from the cosmic shear signal. We then 
establish an estimator to extract the galaxy density-density-intrinsic ellipticity correlation (ggl) from the observed galaxy 
ellipticity-density-density measurement for a photo-z galaxy sample. Section 4 describes the residual sources of error to the 
Gil self-calibration technique, and we present the necessary relations to quantify these errors. Finally, we summarise the 
effectiveness and impact of the Gil self-calibration in Sec. 5. We expand in the Appendix upon the detailed calculation of the 
coefficients in the error calculation found in Sec. 4.1 and provide a list of expected values. 



2 BACKGROUND 

The self-calibration technique proposed by Zhang (2010a) for the 2-point intrinsic alignment correlations and generalized 
by Troxel & Ishak (2012a) for the 3-point correlations makes use of the information already found in a lensing survey, 
including galaxy shape, angular position and photometric redshift, in order to calculate and remove the dominant intrinsic 
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alignment contaminations. In evaluating the performance of the self-calibration technique, we will consider as an example 
survey parameters to match an LSST-like weak lensing survey (LSST Science Collaborations and LSST Project 2009), but of 
course the calculations are applicable to all current and planned weak lensing surveys (e.g. CFHTLS, DES, EUCLID, HSC, 
HST, JWST, LSST, Pan-STARRS, and WFIRST). We divide galaxies into photo-z bins according to photo-z z p , where the 
i-th photo-z bin has a range z\ — Az^/2 ^ z p ^ Zi + Azi/2 for mean photo-z z%. In this notation, i < j implies that Zi < Zj. 
The galaxy redshift distribution over the i-th redshift bin is n p (z p ) and rii(z) as a function of photo-z and true redshift, 
respectively. These are then related by a photo-z probability distribution function (PDF) p(z\z p ). 

Our calculations assume that the survey will acheive a half sky coverage (f s k y = 0.5) and a galaxy surface density of 40 
arcminute -2 with redshift density distribution of 

n(z) = JT ("§ ) ex P(-^)' W 

with zq = 0.5. The ellipticity shape noise is described by 7 rms = 0.18 + 0.042z (Zhang 2010a) and the photo-z error by a 
Gaussian PDF of the form 

^/2^f P V 2a 2 2 

with a, = 0.05(1 + z). We construct photo-z redshift bins with width Az = 0.2, centred at 2, = 0.2(i + 1) (i = 1, • • • , 9). 
Redshifts below z p = 0.3 are excluded, not because of poor performance in the Gil self-calibration technique, but rather due 
to the weaker lensing signal at low redshifts, which artificially inflates the fractional errors we evaluate in Sec. 4. These errors 
are measured with respect to the lensing signal, and thus lower photo-z bins are not useful in evaluating the true performance 
of the Gil self-calibration. 

The intrinsic alignment self-calibration techniques (Gil and CGI) both rely upon two basic observables measured in a 
weak lensing survey: galaxy surface density and galaxy shape. The galaxy surface density, <5 E , is a function of the 3D galaxy 
distribution 8 g in a given photo-z bin. The galaxy shape is expressed in terms of ellipticity, which measures the cosmic shear 7. 
However, the intrinsic ellipticities of galaxies contaminate the cosmic shear. This intrinsic ellipticity has a random component, 
which is simple to correct for and which we include as part of the shot noise in the error estimations of Sec. 4. A second 
component to this intrinsic ellipticity is due to the correlated intrinsic alignment of galaxies and was introduced in Sec. 1. 
The measured shear can be labelled as 7 s = 7 + 7 7 , where denotes the correlated intrinsic ellipticity due to intrinsic galaxy 
alignment. We are concerned only with the weak limit of gravitational lensing, so we can work with the lensing convergence 
k instead. Thus from the measured 7 s , we can obtain the convergence k s = k + k 1 . 

We assume a standard, fiat ACDM universe in our calculations. The convergence k of a source galaxy at comoving 
distance \G and direction is then related in the Born approximation to the matter density 5 by the lensing kernel Wl(z' , z) 

k{0) = f XG S( XL , 8)W L ( X L,XG)d X L. (3) 
Jo 

The 3D matter bispectrum is defined from k as 

(k{£i)k(£ 2 )R(£ 3 )} = (2Tr) 2 S D (£ 1 +£ 2 +i3)B 6 (£ 1 ,£2,e3), (4) 

where the ensemble average is denoted by (• ■ • ) and 5 D {£) is the Dirac delta function. S D (£1 +£2 +£3) enforces the condition 
that the three vectors form a triangle in Fourier space. Under the Limber approximation, we express the 2D angular cross- 
correlation bispectrum as 

B^(£;z[,z p ,z p ) = f X W °^ ; *' Xa -* 3) iW*;xW, (5) 
Jo X 

where when a = f3 = 7 = G, for example, Bggg {k; x') is the 3D matter bispectrum shown in Eq. 4. However, more generally 
a, f3, 7 G G, I, g, and the intrinsic alignment (I) and galaxy (g) bispectra are calculated as described below. The redshift is 
simply related through the Hubble parameter, H(z), to Xi an d we can write the weighting function in terms of redshift as 

^ 7 ( z (x)^r(xi),^ P (X2),^(X3)) = Wt(z,z P )Wf(z,z[)W^(z,z P ), (6) 

TO I P\ I TXT I I \ P 



Wi(z,z P ) = / W L (z',z)n P {z P )dz, (7) 
Jo 

W?(z,z p ) = Wf(z,z p ) = n p (z p ). (8) 

Wp(z, z p ) is simply the weighted lensing kernel for redshift bin i, while Wi(z,z p ) = W®(z,z p ) is the normalized 
galaxy photo-z distribution. The required lensing, intrinsic alignment, and galaxy bispectra are calculated following Troxel 
& Ishak (2012a) and Troxel & Ishak (2012b). A deterministic galaxy bias is assumed for the galaxy bispectra, and the 
galaxy intrinsic alignment signal is calculated using the fiducial parameterized model of Schneider & Bridle (2010) (SB10), 
which is based on the halo model prescription. The SB 10 model reduces by design to the linear alignment model (Hirata & 
Seljak 2004) at large scale, but aims for a more motivated modelling of intrinsic alignment at small scales. 
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The flat-sky bispectrum is related to the all-sky bispectrum through the Wigner-3j symbol, where 

o o 3 ) ^^s^+n^fc,*,. m 

Following the fitting formula of Scoccimarro & Couchman (2001) with coefficients F 2 B (ki, k 2 ) described in Section 2.4.3 of 
Takada & Jain (2004), we compute the 3D bispectrum due to nonlinear gravitational clustering 

B s (ki,k 2 , k 3 ; X ) = 2Ff(k 1 ,k 2 )Ps(k 1 ; X )Ps(k2; x) + 2 perm. (10) 
We make a direct expansion of this method to approximate the 3D intrinsic alignment bispectra, using the intrinsic alignment 
power spectra instead of the nonlinear matter power spectrum, where 

B«,i(*i,fc,fc3;x) = 2Ff(k 1 ,k 2 )P s¥ (k 1 -x)P6(k 2 ; X )+'2Ff(k 2 ,k 3 )P s (k 2 - X )P s¥ (k 3 -x) (11) 

+2Ff(k 3 , ki)P s¥ (fc 3 ; x)Ps-y (ki;x) 
B s¥¥ (kuk 2 ,k 3 ;x) = 2Ff{k 1 ,k 2 )P ¥ {kr,x)P S¥ {k 2 -x) + 2Ff(k 2 ,k 3 )P s¥ {k 2 - X )P s¥ {k 3 -,x) (12) 
+2F 2 cff (fc 3 , fei)P^ (As; x)Ps-V (ki;x) 
B ¥ (ki,k2,k 3 ;x) = 2F 2 cff (fc 1 ,fc 2 )P^(fci;x)^(fe;x) + 2 perm. (13) 
This treatment provides reasonable results for the intrinsic alignment bispectra. 



3 THE Gil SELF-CALIBRATION TECHNIQUE 

A lensing survey will capture the information necessary for several sets of correlations between the 5 g and k s which can be 
constructed for galaxy triplets. Like the GGI self-calibration technique, the Gil self-calibration requires only three of these 
observed correlations. The first is the angular cross-correlation bispectrum between galaxy ellipticities (k s ) in the i-th, j-th 
and fc-th redshift bins 

B$(*L,*a,/ S ) = B^ GG (h,£ 2 ,£ 3 ) + Blf k G {lx,i 2 ,t 3 ) + (2 perm.) + B 1 ^ {£i,£ 2 ,£ 3 ) + (2 perm.) + B 1 ^ (h,£ 2 ,£ 3 ). (14) 

Unless catastrophic photo-z errors overwhelm the data, we can safely neglect the III cross-correlation, which is negligible 
for thick photo-z bins, by selecting galaxy triplets where i < j < k. Though the impact of the Gil correlation under these 
conditions was neglected in Troxel & Ishak (2012a), it can still impact the bispectrum for adjacent bins and should be taken 
into account if the goal is to make a precise measurement of the total intrinsic alignment contamination. For bin combinations 
which are not adjacent, this Gil term can again be safely neglected and the GGI self-calibration is sufficient. Under this 
requirement, we still find that B\j G 3> B 1 ^ , B GI k ' due to the lensing geometry. We then have for i < j < k, 

B$(ii,*2,is) « B? GG (£ 1 ,£ 2 ,£ 3 ) + B^ k G {£ u £ 2 ,£ 3 ) + Bgf(t lt t2 t ta). (15) 

The contribution from the Gil bispectrum B^j G (£i,£ 2 ,£ 3 ) (i < j < k) is typically very small and in some cases effectively zero. 
The Gil self-calibration technique, which seeks to calculate and remove the Gil cross-correlation which survives in adjacent 
photo-z bins, then acts as a correction to the GGI self-calibration, which is the only intrinsic alignment contamination 
that survives for non-adjacent photo-z bin combinations. We will denote the Gil bispectrum B\j G in order to preserve the 
association of each quantity G or / to its redshift bin. 

The second observable is the angular cross-correlation bispectrum between convergence (k s ) in the i-th redshift bin and 
the galaxy density (5 s ) in the j-th and fc-th redshift bins. The self-calibration requires the case where i — j = k 

B%(£i,£2,£ 3 ) = B GGa (£i, £ 2 ,£ 3 ) + S G / fl (£!,£ 2 ,£ 3 ) + B%?(£ x ,£ 2 ,£ 3 ). (16) 

This correlation contributes further information about the intrinsic alignment of galaxies. 

The final observable is measured in the angular cross-correlation bispectrum between galaxy density (<5 S ) in the i-th, j-th 
and fc-th redshift bins for i — j — k, giving 

B%(tui2,£ 3 ) = B™ 3 (£ 1 ,£ 2 ,£ 3 ). (17) 

We also require for the Gil self-calibration results from both the GI cross-correlation and intrinsic alignment auto-correlation 
self-calibration techniques (Zhang 2010a,b). We have neglected thus far the contribution of magnification bias to these 
measurements. This was discussed and justified for the 3-point measurements in Sec. 4.3 of Troxel & Ishak (2012a) and 
was shown to be negligible. There is also a non-Gaussian contribution to the observed bispectra, which is also neglected as 
discussed in Sec. 4.4 of Troxel & Ishak (2012a). 

Our Gil self-calibration technique calculates and removes the Gil contamination in Eq. 15 through the measurements 
from Eqs. 16 & 17, which are both available in the same lensing survey. The fractional contamination of the correlated intrinsic 
alignment to the lensing signal is expressed as 

f i ( p e n \ - B " G (£i,£ 2 ,£ 3 ) 

f^ iM3 ) = -^-—. (18) 
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For the self-calibration to work, the contamination fL^[li,l2,ls) must be sufficiently large as to contribute a detectable 
B^ s (£i,£ a ,ta) at the corresponding I bins in 5^(4,4,4). We denote this threshold f% k resh . When f? jk ^ f!!% esh , the 
Gil self-calibration can be applied to reduce the Gil contamination. The residual error after the Gil self-calibration will 
be expressed as a residual fractional error on the lensing measurement. We differentiate Afij k as statistical error and 5/y/t 
as systematic error, as in the discussion of the GGI self-calibration. The performance of the Gil self-calibration will then 
be quantified by the parameters f*j£ , A/y/t and Sfijki which are discussed and calculated in Sec. 4. For parallelism, we 
preserve the notation from Troxel & Ishak (2012a), though it is essential to note that the quantities are often not the same 
in each technique. 



3.1 Relating B// fc G and B^ g 

In order to self-calibrate the Gil bispectrum, we must first determine the relationship between B//j? and B^f . We follow 
the same general approach as for the GGI bispectrum, relating the two by use of a deterministic galaxy bias b 9tk (Fry & 
Gaztanaga 1993) such that the smoothed galaxy density is a function of matter density expressed as 

5 g (x; X ) = b 3 Ax)S m (x; X ) + ^^8 2 m (x; X ) + 0(S 3 ). (19) 

The linear galaxy bias (6 s .i) was used by Zhang (2010a) for the 2-point correlations. In addition, b 3t2 represents the first order 
non-linear contribution. b 9:2 is typically expected to be negative and ^ b g .i (Cooray & Sheth 2002), and it is insufficient to 
model the bias as simply scale dependent as in the 2-point case (Jeong & Komatsu 2009). Following Troxel & Ishak (2012a), 
we use this expression for the galaxy density to relate B^^i to B^ lB I . We neglect the portion of the bispectrum due to 
primordial non-Gaussianity and the trispectrum term, which contains further information about the non-Gaussianity. This 
was justified and discussed further in Sec. 4.4 of Troxel & Ishak (2012a), and results in the relationship 

B gl i 1 i(k 1 ,k 2 ,k 3 ;x) = b g,i(x) B s-,' (kl, k 2 , k 3 ; X ) 

+ b g ,i(x)b a Ax) [p 7777 ( fci;x )p i7/ ( fc2;x ) + p a7 ,( fe;x )p 577 ( fe;x ) + p 7j77 ( fci;x )p ay ( fe;x )] . (20) 

If the galaxy bias changes slowly over the i-th redshift bin with median comoving distance Xs, we can write to a good 
approximation b\ — b 3t k(xi)- We can further approximate B(ki,k 2 ,ky,x) ~ B(ki, k 2 , k$; x») and P{k;x) ~ P(k; X i) in the 
limit where the comoving distance distribution of galaxies in the i-th redshift bin is narrow, and substituting Eq. 20 into 
Limber's approximation for B 1 ^ 9 and comparing to B^^ , we have the approximations 

B${li,l*,h) « ^ (biB^i (k u k 2 ,k 3 ; X i) 



b\bl 



and 



[Pyi~/i {kv,Xi) p s^ {k 2 ; Xi ) + P Sl i (k 2 ; Xi) p s^ ( k r, x0 + fyy {kv, Xi) p s-y' ( k a; xO] ) > (21) 



B{£ {I x , i 2 , ta ) « B^i 7 / (fei , k 2 , k :i ■ X i) ^% ■ (22) 

Xi 



where W ijk = J °° fi(x) fj{x)W k {x)dx an d TL%h = f£° fi{x) d X- From Eqs. 22 & 21, we can now write 



B ijk («1,*2,«3) ~ , in ±1 iH (ll, 12,13) —J 

Ojiiiii Xi 

x [P^ii^^x^Pssik^Xi) + Pss(k 2 ;xi)P^yI(k3;x^) + Psy'(ki;Xi)Ps~,'(k:i;Xi)] ■ (23) 



In order to express the 3D power spectra in Eq. 23 as 2D spectra, we will use the approximations C^f(£) ~ P s i (k; Xi)~, 



bi n, 



and C"{€) « P 7 7 7 /(fc;xi)^¥-- where Uu = J °° ff(x)dx- Equation 23 is then 

BiffdMs) ~ ^B»°(l lM )-ti^ 

x [cli{h)Cl?{l 2 ) + KC*f(t a )Ctf(ta) + CU{£i)Cl?{h)] . (24) 

This expression now relates B\j k to B^ g , allowing us to identify the Gil bispectrum in the observable B^ 2 \ Like the 
similar expression for the GGI self-calibration, it is necessary to have information about the gl and II power spectra obtained 
from the 2-point self-calibration techniques (Zhang 2010a,b). 
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(B«? S |s + B^ 3 \s)/(Bf 1 T 3 + B^l 9 ) for equilateral triangles (I = t x = l 2 = £ 3 ) over three redshift 



bins spanning the survey range. The expected values of Qgii arc similar to those found for the GGI self-calibration, though with 
more ^-dependence. This is primarily due to the greater complexity found in Qgii, as the sum of two components instead of one. The 
suppression is dependent on the redshift bin chosen, increasing with redshift due to increased photo-z error at higher redshift. Generally, 
Q a 0.4, and the significant deviation from unity ensures that the estimator B^ 9 is valid for lensing surveys of interest. 



3.2 B. 



iig 



estimator 



We can measure B\l 9 through the information contained within the observable = B\^ 9 + B^ 9 + B 



GGg 



To measure 



it directly, we must first remove the contamination of the geometry dependent BfJ 9 + B^ 99 , which is not useful for the 
Gil self-calibration. Lensing geometry simply requires eliminating those triplets of galaxies where the redshift of the galaxy 
used to measure the ellipticity is lower than those used to measure galaxy number density in the case of a spectroscopic 
galaxy sample. In this way, those triplets remaining have no contamination from lensing and measure only B\l 9 . However, 
for a photo-z galaxy sample, the typically large photo-z error prevents us from directly employing this method. Photo-z error 
causes a true redshift distribution of width ^ 2o\p = 0.1(1 + 2) even for a photo-z bin with width Az — > 0. In practice, photo-z 
bin widths are typically ^ 0.2. It is possible with such photo-z errors for galaxy triplets in the i-th redshift bin to provide a 
measureable lensing contribution even when i < j, k, except for the special cases where the redshift or both the photo-z error 
and bin size are limited to sufficiently low values. A general photo-z galaxy sample then requires a more careful approach 
when separating B ! J 9 from B% 9 + B<f 99 . 

In order to construct such an estimator, we first consider the orientation dependence of the two components. We define a 
redshift for each galaxy in the triplet: zg, Zq', Zi, or zp for the galaxies used in the lensing/intrinsic alignment measurement 
and Zg for the galaxy used in the number density measurement. The gll correlation is independent of the relative position of 
the three galaxies. For example, the correlations with zi < Zji < z g , z g < Zi < Zv or zp < z g < zi are statistically identical 
when the sides of the triangle are fixed. However, the GGg and Gig correlations do depend on the relative position of the 
three galaxies. Due to the lensing geometry dependence, the correlation with zg, zq' ii < z g is statistically smaller than other 



orientations. 

This dependence provides two observables from . 



The first is B^ , where all triplets are weighted equally. The second 
which counts only those triplets with zg,zg> < z g . This weighting is denoted by the subscript 'S'. This ordering is 
necessarily different from that employed in the GGI self-calibration. From our previous discussion, we then have B\l 9 = B ! J 9 \s 



is B 



(2), 



and B^ 9 + B^ 9 > \B<^ 9 + B% 3 ]\s = B^ 9 \ s + B^ 



We now define the suppression ratio 







B^ 9 \ s (£i,i2,£ 3 ) + B?J 9 \ s (£i,l2,l 3 ) 



B 



GGg 



(£i,e2,e 3 ) + B^ 9 (£ 1 ,£ 2 ,£ 3 



(25) 



where we have explicitly included the ^-dependence which had been previously neglected. This ratio describes the suppression 
of the signal due to the weighting of triplets described for subscript 'S'. By definition < Qgii < 1, with Qgii — if the 
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photo-z is perfectly accurate and Qgii = 1 if the photo-z has no correlation to the true redshift. Qgii is calculated using the 
galaxy redshift distribution, which is discussed in Sec. 3.3. 

The definition of Qgii(£i, £2, £3) possesses an explicit dependency on the intrinsic alignment, which in principle would 
make the estimator dependent on intrinsic alignment model. However, as we discuss in Sec. 3.3, we find that for an LSST-like 
survey, B GGa is large enough compared to B G ? g that we can approximate Qgii(£i, £2, £3) ~ B GGs \s /B GGg to very high 
accuracy. We can thus treat Qgii(£i, £2, £3) as independent of intrinsic alignment model for sufficiently deep surveys. 

We now define the estimator for B^ 3 , which we denote B^f , in terms of Qgii(£i, £2, £3) and the two observables 

B%(£ 1 ,£ 2 ,£3) = B I J/(£ 1 ,£ 2 ,£3) + BfJ 3 (£ 1 ,£ 2 ,£ 3 ) + B^ 3 (£ 1 ,£ 2 ,£ 3 ), 
B^siluM) = B I J/(£ 1 ,£ 2 ,£3) + BfJ 3 \s(£i,£2,£3) + Bfif !3 \s(£i,£2,£3). (26) 



This estimator is 



VH,,. v. ^- B«i\s(li,t2,h) -Qgii(£i,£2,£3)B%( 



Bm (£i,£2,£3) = : , 9 \ ■ i - 27 > 

1 — Qgii(£i,£2,£3) 

As expected, when Qgii = this gives B\l 3 — Bj?) \s as for a spectroscopic galaxy sample with no photo-z error. However, 
Qgii must not approach unity, where B\l 3 is singular. For the LSST-like survey described in Sec. 2, we calculate Qgii for 
various redshift bins following the procedure described in Sec. 3.3. This result is given in Fig. 2 for equilateral triangles, where 
we find Qgii ~ 0.4 and in general that Qgii should deviate significantly from unity. The estimator B\l 3 is thus expected to 
be applicable in any typical lensing survey. 



3.3 Evaluating Qgii(£i, 4, 4) 

The ratio Qgii in Eq. 25 must first be evaluated in order to employ the estimator B ! J 3 . This mirrors the derivation in Sec. 
3.3 of Troxel & Ishak (2012a), which we will summarize here as it applies to Qgii- We begin from the real space angular 
correlation functions w GG 3 (0i, 6 2 , 83; Zq, Zq, , z P ) and w GIg (#1, 6 2 , 83; Zq, zf , z P ) . Taking the average correlation, we express 
this in terms of the ensemble average and calculate the associated bispectra. For the correlations with galaxy triplets weighted 
by the subscript 'S', we define the statistical weighting functions 

, > 3/, dz p G J. dzf J. dz p g J °° dz G W L (z L , z G )N p S(z£, zf, z P ) 

n[ZL > Zl ' Za) J i dzSJ i dzfJ i dzar^GW L ( Z L,ZG)Nf {2 *> 

3/. rizg /. dz P J. dz p , J °° dz G J °° dz G 'W L {z L , z g )W l (zl, z G ')N' P S(z P , z p ,,z p ) 

V [ZL ' Z9 ' } I, dz P J, dzf J. dz P /~ dz G J" dZG l W L {z L ,ZG)W L {z L ,Z GI )N' P 1 > 

where 

N P = p(z G \zG)p(z I \zf)p(z g \z P )n P (zG)n P (zf)n P (z P ) (30) 
N'f = p(zG\zc)p(zG'\zG')p(z g \z P )n P (zo)n P (za')n P (z P ), (31) 

and S(zq, zf , z P ) — 1 {S(z G , z g i , z P ) = 1) if z G zf < z P [z g ,Zqi < z P ) and is zero otherwise. Since S(- ■ ■ ) allow only 1/3 
of the integral to survive, r], r\ are normalised by a factor 3 in order to remove the suppression due to the selection function 
and measure only that due to the lensing geometry. 

We can now express directly the bispectra necessary to compute Qgii 

B% 3 (£ U £2,£ 3 ) = rB Gl3 (k 1 ,k 2 ,k3;x) m{x)f : {x) d X (32) 
Jo X 

B GG3 (£ 1 ,£ 2 ,£ 3 ) = I" B GG °(k u k 2 Mlx) Whx) l i{x) dx (33) 
Jo X 

and 



B G J 3 \ S (£ 1 ,£ 2 ,£ 3 ) = r B Gl3 (k 1 ,k 2 ,k 3 ; x ) Wi{X }l\ {X \ {x, x(*i) = X, x(z B ) = x)d X (34) 
Jo X \ z ) 

B GG3 \ s {h,£2,£3) = r B GG3 (k u k 2 ,k 3 ; x ) ^M^ v'U, xM = X) = x)*X- (35) 
Jo X \ z ) 

The ratio Qgii is now defined through Eqs. 33 & 35. For the deep survey we are considering, we find that we can 
approximate to very high accuracy Qgii ~ f]' \, where rj' i is the mean value of rj across the i-th redshift bin. This is because 
for a deep survey, B GGg is much larger than B G ' 3 , making the contribution from 77 negligible. In the limit where photo-z 
error dominates, up 3> Az, and so 77, 77' — > 1. In this limit, the estimator B\l? becomes singular and B\l 3 can no longer be 
differentiated from B GGg + B G ' 3 . In the opposite limit, where a v -C Az, r\,r\ — > 0, where our estimator mirrors the extraction 
method for B\l 3 in spectroscopic galaxy samples. 
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4 PERFORMANCE OF THE Gil SELF-CALIBRATION 

We will quantify the performance of the Gil self-calibration technique using the survey parameters described in Sec. 2. 
This includes both a statistical measurement error introduced through the estimator Bf u 3 , which propagates into the final 
measurement of -8//?? and thus of B GGG , and a systematic error due to inaccuracy in the scaling relation between BU g 
and B\l? . We will also summarize other possible sources of error which might impact the performance of the self-calibration, 
though a more detailed account can be found in Troxel & Ishak (2012a). 



4.1 The estimator B" 3 

The estimator B\l 3 introduces a statistical error into the measurement of B\* 3 and thus B GGG . In order to quantify the 
accuracy of the estimator B^ u 3 (£i, £2, £3), we consider the contribution of measurement errors such as shot and shape noise 
m B^(ei,e 2 ,£ 3 ) which propagate into our measurement of B^ 3 (£\, £2, £3) through the estimator. For a given redshift bin 
we calculate the rms error, working in a pixel space with Np sufficiently fine and uniform pixels of photo-z z& and angular 
position 9 a . The measured overdensity is given by 5 a + 5^ , and the measured 'shear' is given by K a + K ! a + , where 'N' 
represents the measurement noise. From Eq. 26, we construct the pixel space angular bispectra 

B i2) (£ 1: £ 2 ,£ 3 ) = Np 3 J2[S a + S^][np + 4 + + «4 + K?]e*p[t(/i • a +£ 2 Op +l 3 ■ fty)]. 

B (2) \ s (£i,£ 2 ,£ 3 ) = Np 3 Y^i s c'+S^][^ + 4 + 4][^+4 + ^] e MK£i-Ooc+£2-0p+£3-e^)]S a ^. (36) 

ck/3 7 

S a p-y = 1 when 2q > zp , z^ and is zero otherwise. Thus in the limit Np 3> 1, X] a /3 7 Saj87 = Np/3 and the average S a p~, = 1/3. 
From our definition of the estimator in Eq. 27, we can construct the difference 

B&'-Bg' = — -1 riVp 3 ^exp[i(£ 1 .0 Q +£ 2 .0^+£ 3 -^)][(3^ 7 -Q G/I ) 

1 1 1 



X(S a + 5 a )(Kp + Kp + K^)(K 7 + K 7 + K y ) — (1 — Qai^&aKpKy 

= 7- ^ TNp 3 V exp[i(£i ■ a + £2 ■ Op + £3 ■ 7 )](3S a(37 - Q 3 ) 

(1 - Qgii) "V 

X [((5c + 5a)((Kf3 + hip + Kp)(H-y + ) + (up + )« 7 )]- (37) 

We have used here that <S a /3 7 = 1/3 and that the gll correlation doesn't depend on the relative position of the galaxy 
triplets. We can now write the rms error in a similar way to Eq. 41 of Troxel & Ishak (2012a) and simplify the resulting 
256 6-point correlations which arise from expanding Eq. 37. We apply Wick's theorem, which allows us to express each 
6-point correlation as 15 products of three 2-point correlations. This results in 3840 products, most of which are zero. For 
example, any correlation between signal and noise or dissimilar noise terms vanish. Further, due to the angular dependence 
of the correlations ((A a Bb) = wab(Q<i — 9b)), only those correlations with (A a Bb) where a £ a,P,j and b £ A, /i, v are 
non-vanishing. This leaves 124 surviving products in the rms error expression 

(AB'jA 2 = - -j Np 6 V V exp [i{t x O a +£2-O +£3- fty)] cxp [i(£x ■ X + £ 2 ■ 0„ + £ a ■ ft,)] (3S Q/37 - Q 3 ) 

X (3Sa p „ - Q3) I + (k^«m )) [((^q^a) + (6a 8\)) (((«7 + I-y)(Kv + Iv)) + (k*k*)) 

+ (S c ,(k 1/ + Z„))((k 7 + / 7 )5a)J + ({^pK v ) + (KpK^)) U(5 a 5x) + (8a8x)) + + ^)) + s)) 

+ + /m)><(^7 + I-y)Sx)] + ««7«m> + <«» [«W + <W» «(«/> + W> + (Jp/C„» 

+ + I*)){{K P + Ip)Sx)] + + ««?» [«<Wa> + <<«>) «(«/> + /|))^> + (38) 

+(5 Q (^, I + I m ))((k/3 + ^)<5a)] + ((5«5a) + (Sa^x)) Unpl,*) (I-yiu,) + (npI v )(IyKn) + (ttyln)(lpn v ) 
+(^ 7 J,y)(J ( 3K M )j + (Kp5\} ^(5 a K M )((K 7 + 7 7 )/ y ) + (5 Q (fi: M + 7 M ))(7 7 /t I/ ) + (<5 a K I /)((/t 7 + / 7 )/ M ) 

+ (5 q (k„ + J I/ ))(/ 7 K M )j + (k 7 (5a) [(<$ a KVi)((/«/9 + 7/3)^) + (<5aKi/)((«^ + Ifi)In) + (<5a(«i' + 

+ (5q(^, i + J M ))(//3K I /)j + (//3<5a) [(5aKf4)( K 7^) + (^a«i-) + (^7^) [('JaKfj) (Kfil v ) + {S a K, v ) {hZpI ^ j. 

Noises only correlate at zero lag ((S^Sx) oc S aX , {k£ 1$) oc 5 7I/ ), and the correlations not involving k depend only on 
separation, not on relative orientation of the galaxy pairs along the line-of-sight. However, correlations like (n6) and (kk 1 ) 
are dependent on the relative orientation and must be treated with care when evaluating Eq. 39. In order to quantify this 
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le-05 1 1 1 1 1 1 1 — 1 — 1 — 1 1 

100 1000 

/ 

Figure 3. The residual statistical uncertainty A/|?2 for equilateral triangles [t = t\ — = I3) in the Bff^ measurement and threshold 
of intrinsic alignment contamination /*JjT es at which the Gil self-calibration technique can calculate and remove the intrinsic alignment 
contamination at S/N=l arc plotted for a variety of redshift bin combinations. At large £, we see the effects of shot noise beginning to 
take over. Generally, is less than the minimum survey error expected for such a survey, and is thus negligible. We expect this 

result to hold for non-equilateral triangles as well, but the use of the Gil self-calibration is limited by our understanding of non-Gaussian 
effects for very elongated triangle shapes, as discussed in Troxel & Ishak (2012a), and we leave discussion of its applicability for these 
very elongated triangle shapes to a future work. 



orientation dependence, we apply Qi g = Q2 (Zhang 2010a; Troxel & Ishak 2012a) and Qig such that 

/ c \ 1 / Socu Sua \ 1 e \ 

^-K(^fcy + fe)^- (39) 

The suppression ratio Qig ~ Qig is denned in an identical way to Qi g , but for the C IG power spectrum. We can now 
evaluate Eq. 39 analytically, converting to its Fourier representation 

(AB^y = 2{ (C GG + C u ) C Gs C Gs Kx + C GG [ (C GG + 2C n ) (c gg + C ggI V) + 2C Ig C Ig 
+C Gg C lg K 2 + C lG (c gs K 3 + C* ggN A' 4 ) + C* GGN (c gg A 5 + C ggN A 6 ) ] 

+C IG r^IG / C gg^ 7 + C ^ Kg \ + C G SC G SK9 + C Gs C IS Kl0 (40) 



+C GGN f C SS Kll + C SSN K 



12 



+C GGN J C GGN ^ c ,gg Ci + C ggN ei \ + c^S^ + qH (c^ Kb + C ^ N K 6 

+C Gg C* Gg A 13 + C Gg C Ig A 14 ] } 

The details of this calculation and the coefficients K\ — A14 and ax, cx, e\ are included in the appendix. 
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The final rms error 

AB 'iu evaluated for a given triangle with bin width A£ is then given by 

K»-r - ,,,,,, M i aw j ^ tci ')^ 

+C GG [ (C GG + 2C n ) (c gg + C sst V) + 2C Ig C Ig 

+C Gg C Ig A'2 + C IG (c gg A" 3 + C ggN A 4 ) + C GGN (c gg A 5 + C ggN A- 6 ) ] 

+C IG [C IG (c gg A 7 + C ggN K 8 ) + C Gg C Gg A 9 + C Gg C Ig K 10 (41) 
+C GGN (c^Kn + C ggN A 12 ) ] 

+C GGN [C GGN (c gg Cl + C ggN ei ) + C Ig C Ig A 5 + C n (c gg A 5 + C ggN A- 6 ) 
+C Gg C Gg A' 1 3 + C Gg C Ig A 14 ] } 

Cfi' N = and C GG ' N = yrms/^i, where Hi is the average number density of galaxies in the i-th redshift bin. The statistical 
error AB^f differs in several key ways from the statistical error AB^f in the GGI self-calibration, both in the complexity of 
its calculation due to the additional k term and the resulting components, which include terms like C IG . It is more sensitive 
to the intrinsic alignment contamination in the limit than in the Igg and 2-point cases, but is still safe from strong dependence 
in the limit where C GG is dominant. 

The errors AB^f and AC-f (Zhang 2010a) propagate into the measurement of B\l G through Eq. 24. To find the fractional 
error Af^ this induces in the lensing bispectrum, we simply scale AB-j k G by the factor /yj. such that A/^ = A-B^-jF. 
This error is equal to fij k eeh , the minimum intrinsic alignment f?j k which can be detected through the self-calibration with 
S/N=l or AB\l? = B"?. Thus both the residual statistical error in the measurement of C GGG and the lower limit at which 
the intrinsic alignment can be calculated and removed with the self-calibration is represented by fijk esh = Af^ k . Thus the Gil 
self-calibration technique can turn a systematic contamination /y fe of the lensing signal into a statistical error Af^ < flj k , 
which is insensitive to the original intrinsic alignment contamination. 

We show in Fig. 3 Af^ k for the survey parameters described in Sec. 2. Compared to the minimum survey error as found 
in Troxel & Ishak (2012a), A/y^ is negligible for most scales, only becoming comparable for some photo-z bin combinations 
at very large scale. 



4.2 The accuracy of the B^-B'J/ relation 

In addition to the statistical error introduced through the estimator B^f , there is a systematic error which is introduced by 
Eq. 24, which relates the intrinsic alignment contamination B\j h G in the lensing bispectrum to other survey observables. The 
accuracy of Eq. 24 is quantified by 

W, jk B I l ^(£ 1 ,£ 2 ,£ 3 ) 



1 !.)'>' 



b\ W ijk 



(&i) 2 n iM B I ij f(£i,£ 2 ,£ 3 ) OA) 2 ujuILu B*jf(ii,t a ,£a) 
Cl?{h)C GG {h) + C GG (£ 2 )Ci 3 (£ 3 ) + -^-Cli^ClHh) 
This induces a residual systematic error in the lensing measurement of 

cr sys rl 

°Jijk — ^ijk Jijk- 



(42) 



(43) 

e ijk ^ S evaluated numerically and shown in Fig. 4 for equilateral triangles and a variety of neighboring photo-z bin choices. 
We neglect totally non-adjacent photo-z bin combinations, as the Gil signal is truly negligible in such cases. The accuracy of 
Eq. 24 reflects the differences between B\j k G and B(j% G . The accuracy of Eq. 24 depends primarily on mean redshift of the 
photo-z bins chosen, becoming more accurate at higher redshift where the intrinsic alignment is less strong. By comparison, 
the accuracy of the equivalent scaling relation for the GGI self-calibration shows strong tendencies toward the effects of the 
lensing kernel, since the lensing contribution to the signal is much stronger. Equation 24 is accurate to within ~ 30% for all 
scales and photo-z bin combinations. However, for most scales and photo-z bin choices, it is accurate to within 10%. This 
is comparable to the accuracy found in the GGI self-calibration, and we expect the Gil self-calibration to be capable of 
reducing the Gil intrinsic alignment contamination by a factor of 10 or so in general. These results are insensitive to the 
original intrinsic alignment contamination, such that for any fij k esh < flj k < 1, the GGI self-calibration will reduce the GGI 
contamination down to survey limits or by a factor of 10 or so, whichever is less, for all but a few photo-z bin choices at large 
scale. The actual impact of the systematic on the measurement of B GGG is significantly less, however, as it is scaled by f'j k , 
which is typically expected to be on the order of 10-20%. This would lead to an actual systematic error on the percent level. 
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Figure 4. The inaccuracy of the relationship between B^ff and the observable B 1 ^ 9 is quantified in Eq. 42 by ejjjj- The dominant 
systematic error in the measurement of B^GG through the Gil self-calibration technique is due to this inaccuracy. For equilateral 
triangles (I = l\ =1-2 = ^3)1 ejjjf is plotted for several sets of adjacent redshift bins, where the stronger dependence of the lensing kernel 
on redshift causes a significantly higher inaccuracy. We neglect totally non-adjacent photo-z bin combinations, as the Gil signal is truly 
negligible in such cases. Equation 24 is accurate to within ~ 30% for all scales and photo-z bin combinations. However, for most scales 
and photo-z bin choices, it is accurate to within 10%. This is comparable to the accuracy found in the GGI self-calibration, and we expect 
the Gil self-calibration to be capable of reducing the Gil intrinsic alignment contamination by a factor of 10 or so in general. These 
results are insensitive to the original intrinsic alignment contamination, such that for any ffj£ e < /yj. < 1, the GGI self-calibration 
will reduce the GGI contamination down to survey limits or by a factor of 10 or so, whichever is less, for all but a few photo-z bin choices 
at large scale. 



4.3 Other sources of uncertainty 

The Gil self-calibration suffers from many additional potential sources of uncertainty, including uncertainties in galaxy bias 
modelling, limitations in the models used for the intrinsic alignment and bispectrum calculations, magnification bias effects, 
and the effects of non-Gaussianity in the bispectrum. Troxel & Ishak (2012a) provides detailed calculations which compare 
the order at which magnification and non-Gaussian effects, as well as uncertainties in the galaxy bias model, might impact 
the self-calibration for the GGI bispectrum. While not all of these calculations apply directly to the Gil self-calibration, it is 
trivial to replicate them following Troxel & Ishak (2012a). They also address a range of other possible factors which might 
influence the self-calibration, including for example cosmological uncertainties and catastrophic photo-z errors. In all cases, 
they find that the dominant impact on the performance of the self-calibration comes from the statistical and systematic errors 
described in Sees. 4.1 and 4.2. We find that this holds for the Gil self-calibration as well, with the exception of uncertainty 
in the galaxy bias parameters, which becomes comparable to the statistical measurement error in B IIG at large scale. 



4.4 Summary of residual errors 

The performance of the Gil self-calibration can be summarised under three regimes, which are defined by the magnitude 
of the Gil contamination as represented by f'j k . The first is where the gll correlation is too small to detect in B^ 2 ' , with 
fijk ^ fijk eah - If the intrinsic alignment cannot be detected in B^ 2 \ the Gil self-calibration is not applicable. This generally 
means that the Gil contamination is also negligible when compared to e™™, the minimum statistical error in the lensing 
bispectrum, and there is no need to correct for it. This is likely true for all totally non-adjacent photo-z bin choices, where 
the Gil signal is naturally negligible as discussed in Troxel & Ishak (2012c). 

However, the Gil contamination to the lensing bispectrum is likely not negligible if /y fc > ft ] h k esh and it must be 
corrected for. The Gil self-calibration is now able to detect and calculate the Gil cross-correlation. In the second regime, 
where Af^ esh > ^fljk, the statistical error induced by measurement error in the estimator B ! J 9 is dominant. This 

error is generally negligible when compared to e™^™, and so in this regime, the Gil self-calibration should perform at the 
statistical limit of the lensing survey. 

Where Af^ k reah < ^j^fijki the systematic error 5 fijk = e ijkfljk due to the relationship between B\j h G and B\l? in Eq. 
24 is dominant. In the case where e|J^ < e^™///,^, is still dominant. Otherwise the Gil self-calibration can suppress 
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the Gil contamination by a factor of 10 or so for all but a few photo-z bin choices at the largest scales. In this case, other 
complementary techniques could be employed to further reduce the Gil contamination down to the statistical limit for the 
lensing survey. 

For example, one such case has been explored by Zhang (2010b) for the 2-point correlations, but such studies of 
the 3-point intrinsic alignment are left to be done. Zhang, Pen & Bernstein (2010) combines the GI self-calibration with 
a photo-z self-calibration to better protect the GI self-calibration against catastrophic photo-z effects. Both methods are 
possible because the GI and GGI self-calibration uses primarily those correlations in one redshift bin to estimate the intrinsic 
alignment, while Zhang (2010b); Zhang, Pen & Bernstein (2010) use those correlations between redshift bins. Others have also 
used information between redshift bins to calibrate the intrinsic alignment contamination in the 2- and 3-point correlations 
(Okumura & Jing 2009; Kirk, Bridle & Schneider 2010; Joachimi & Schneider 2008, 2009; Shi, Joachimi & Schneider 2010; 
Joachimi & Bridle 2010). Such techniques for the 3-point intrinsic alignment correlations should eventually complement the 
Gil and GGI self-calibration techniques for improved reductions in the contamination by the intrinsic alignment in the cosmic 
shear signal, but further exploration is necessary in order to jointly apply these complementary techniques to realistic survey 
conditions. 



5 CONCLUSION 

The 3-point intrinsic alignment correlations (GGI, Gil, III) are expected to strongly contaminate the galaxy lensing bispectrum 
at up to the 20% level. Troxel & Ishak (2012c) showed that while the III correlation can be safely neglected by considering 
only the cross-correlation bispectrum between three different photo-z bins, the Gil cross-correlation remains a contaminant 
for adjacent bins in addition to the GGI cross-correlation and thus must be considered in any self-calibration of adjacent 
photo-z bin combinations. Troxel & Ishak (2012a) first generalized the self-calibration technique to the bispectrum in order 
to calculate and remove the 3-point GGI contamination from the GGG bispectrum. In this work we extend the self-calibration 
to the 3-point Gil cross-correlation in order to measure and remove its remaining contamination from adjacent bin triplets. 

In order to do this, we establish the estimator B^ u 9 to extract the gll correlation from the galaxy ellipticity-ellipticity- 
density measurement for a photo-z galaxy sample. This estimator is expected to be generally applicable to weak lensing surveys 
and reduces to the simple extraction method for spectroscopic galaxy samples at low photo-z error. We then develop a scaling 
relation between the Gil and gll bispectra using the linear and non-linear galaxy bias to relate the galaxy density and cosmic 
shear measurements. We can then calculate and remove the Gil correlation from the GGG bispectrum. While this method is 
in principle applicable to all I and triangle shapes, we maintain the modest restrictions of Troxel & Ishak (2012a) on very 
elongated triangles due to the effects of non-Gaussianity and at very non-linear scales due to limitations in the understanding 
of the galaxy bias model used. Combining the Gil and GGI self-calibration techniques will then allow a complete removal of 
the 3-point intrinsic alignment contamination from the cosmic shear signal. 

The performance of the Gil self-calibration technique is also quantified for a typical weak-lensing survey. The residual 
statistical error due to measurement uncertainty in the estimator B^ 9 is shown to be generally negligible when compared to 
the minimum measurement error in the lensing bispectrum. We consider the systematic error introduced by the relationship 
between B^/jF and B\l 9 , showing that \eijk\ < 0.3 for all photo-z bin choices and all scales, while except for the very largest 
scales, |eyfc| < 0.1. The intrinsic alignment contamination can then be reduce the Gil contamination by a factor of 10 or more 
for all adjacent photo-z bin combinations at t > 300. For larger scales, we find that the Gil contamination can be reduced 
by a factor of 3-5 or more. This will potentially allow the Gil self-calibration to reduce the Gil correlation to the statistical 
limit of the lensing survey, as discussed in Sec. 4.4. 

These results are not strongly sensitive to the original intrinsic alignment contamination, such that for any ft^ esh < 
fijk < 1; the Gil self-calibration can reduce the Gil contamination down to survey limits or by a factor of 10 or so, whichever 
is less, for all but the largest scales. This is comparable to the GGI self-calibration, where for any fif resh < /y < 1, it can 
reduces the GGI contamination down to survey limits or by a factor of 10 or greater for most photo-z bin choices, whichever 
is less. We thus expect the Gil self-calibration to perform well with the GGI self-calibration, and together they promise to be 
an efficient technique to isolate the total 3-point intrinsic alignment signal from the cosmic shear measurement. 
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APPENDIX A: CALCULATION OF COEFFICIENTS IN AB^ 

When evaluating the sum and converting Eq. 39 to Fourier space, the products of the correlations each have a numerical 
coefficient due to the restrictions on redshift ordering. Many, however, are identical due to symmetries, and others are 
numerically equivalent when specific suppression ratios are chosen. The calculation of the unique coefficients K\ — iff 6 and 
oi, ci, ei in Eqs. 4f & 42 are summarised here. The coefficients K\ — K16 axe simply collections of other coefficients, as shown 



A" 
*M 




&3 "I" 03 


J\ , 
1\2 




&2 + t>2 + 2C2 






d>2 + 62 + J2 + <?2 


K 4 




fl2 + 12 + J2 + k 2 


K s 




Oi + bi 


K 6 




Cl + d-L 


K 7 




S3 + h 3 


K e 




«3 + J3 


K 9 




«4 + &4 + C4 + ^4 


Ki 




c-i + d-i+ e :i + f 3 


Ku 




h 2 + 712 + 2 + (j2 + <?2)/2 


Km 




ra + S2 + t2 + U2 


Kl 3 




k 3 + h 


K 14 




h + mi + ?2 + r2 



(Af) 



From these, we group the coefficients by the number of orientation dependent correlations, for example (8k) , are involved in 
their calculation. The first coefficient is trivial, due to products with no noise correlations or orientation dependent correlations. 
We then calculate 



(1-Q 3 ) 2 



a/3 7 X\±v 



For the following terms, each is calculated as in the case of Eq. A2, but limited by some ordering restriction due to a noise 
or orientation dependent correlation. For those terms with no orientation dependent correlations, there exist five unique 
coefficients: 



fll = 1+ 20(1-Q G1I)2 (A3) 

&1 S ^(l-Qon)' (A4) 

Cl = 1 + i71 — X rr~vi ( A5 ) 
8 (f - Qgu) 1 

dl = 1+ sn 75 ^2 ( A6 ) 
8 (f - Qgu) 



The coefficient a\ is achieved through the restrictions S a \, 8p v , and 5-y V ; fei is achieved through Sp^ and 5 7 ^; ci through SaxSpp, 
SaxS-y^, 5pn8~, v , and 8p l ,8 ltJ ,; di through 8 a \8p v and 8 a \8 lu ; and finally ei through 8 a xSpij,8 lv and 8axSp v 8j,j,. 
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For those with one orientation dependent correlations, there exist 24 unique coefficients: 
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hi 
ii 
h 

kl 
h 

mi 

Ul 
<>■> 

1>2 
32 

ri 
si 



1 + 
1 + 
1 + 
1 + 
1 + 
1 + 
1 + 
1 - 
1 + 
1 + 
1 - 
1 + 
1 - 
1 + 
1 - 
1 - 
1 + 
1 + 
1 - 



29 - 18Qi g 


11 - 6Qi g 


40 (1 - Qgii) 2 


8(1 -Qgii) 


17 + 6Qi g 


11 -6Qi g 


40 (1 - Qgii) 2 


8(1 -Qgii) 


49 - 18Qi g 


11 -6Qi g 


80(l-Q G n) 2 


8(1 -Qgii) 


61 - 42Qi G 


1 


80 (1 - Qgii) 2 


1 - Qgii 


7 + 6Qig 


1 


20(1 - Qgii) 2 


1 - Qgii 


19 + 42Q IG 


1 


80 (1 - Qgii) 2 


1 - Qgii 


13 - 6Qig 


1 


20 (1 - Qgii) 2 


1 - Qgii 


3(-7 + 5Qi G ) 


1 


20(1 - Qgii) 2 


1 - Qgii 


3(7 + 4Qi G ) 


1 


40 (1 - Qgii) 2 


1 - Qgii 


3 (2 + 5Qi G ) 


1 


20(1 - Qgii) 2 


1 - Qgii 


3(-11+4QicO 


1 


40(1 -Qgii) 2 


1 - Qgii 


3 (5 + Qi g ) 


11 -6Qi g 


40 (1 - Qgii) 2 


8(1 -Qgii) 


3(-8 + 5Qi g ) 


11 -6Qi g 


40 (1 - Qgii) 2 


8(1 -Qgii) 


3 (3 + 5Qi G ) 


1 


40 (1 - Qgii) 2 


1 - Qgii 


3(-8 + 5Qig) 


1 


40 (1 - Qgii) 2 


1 - Qgii 


9(-2 + Qi g ) 


11 -6Qi g 


20(1 - Qgii) 2 


8(1 -Qgii) 


7 + 6Qig 


1 


20 (1 - Qgii) 2 


1 - Qgii 


3(1 + Qig) 


1 


8(1 -Qgii) 2 


1 - Qgii 


3 (-2 + Qig) 


1 


8(1 -Qgii) 2 


1 - Qgii 



(A8) 
(A9) 
(A10) 
(AH) 
(A12) 
(A13) 
(A14) 
(A15) 
(A16) 
(A17) 
(A18) 
(A19) 
(A20) 
(A21) 
(A22) 
(A23) 
(A24) 
(A25) 
(A26) 



ti = 1 + 



1 



11-2 



1 + 



3(l + 3Qi G ) 
8(1 -Qgii) 2 1 - Qgii 
14 - 13Qig 7 + 2Qi G 



8(1 -Qgii) 



8(1 -Qgii) 



(A27) 
(A28) 



Using the shorthand xif = ("W(! ~ Qab) + S p JQab) /2 for Eq. 39, we can express the coefficient conditions as follows. 
The coefficient ai is achieved through the restriction Xc^' ^ 2 l& achieved through Xa?-j c 2 through X\ 9 p and X\\\ di through 
and x\% e2 through xj£ and x 7 ?; h through xf$ and x^; 92 through xl$ and x'£\ h 2 through S^xXf^, 3 a xx'^, $~fvXp°, 
and S^x 1 ^; 12 through 5 aX x'i3° and <5 a AX 7 °; J 2 through 5 a \xi° and 8 a xx'° n and "Vx^a! k 2 through & a \xl$ and <5 q ax'?*; 
h through Sp^x'oH and S^^xLi'-, m 2 through <Wxa 7 and ^cX^i n 2 through Sp^x'y? and $~/vXpZ; 02 through Sp^xly and 
$hj.xI$; P2 through S^xi 9 ^ and 5^xi%i 92 through 5p v Xx% s ~t"X%, and 5p u xi%; T i through <La<Wx 7 ? and <5 q a5 7M x^; s 2 
through SaxSp^xl^ and 5 a \8 ltl xl$\ h through SaxSp^xLj and S aX S 1 „xi$\ and finally u 2 through 8 a \8p v Xw and SaxS^x'pt- 
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Coeflf. Value Coeff. Value Coeff. Value Coeflf. Value Coeff. Value 
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Table Al. Numerical values for coefficients in the calculation of Eqs. 41 & 42 for typical suppression ratios Qgii ~ 2/5, Qig ~ Qig 
1/2. 



For those with two orientation dependent correlations, there exist 12 unique coefficients: 



"3 




1 - 


3 (-24 - 41Qi g + 50Q? g ) 37 - 13Qi g - 2Q? g 


(A29) 




80(1 -Qgii) 2 20(1 -Qgii) 


63 




1 + 


64 - 27Qi g + 28Q 2 g 72 - 25Qi g - 3Q 2 g 


(A30) 




80(1 -Qgii) 2 40(1 -Qgii) 


C3 




1 - 


3 (-23 - Qi g + 3Qi G + 8Qi g Qi G ) 68 - 16Qi g - Qig + 2Q Tg Qi G 


(A31) 




80(1 -Qgii) 2 40(1 -Qgii) 


d 3 




1 + 


3 (20 - 4Qi g + 3Qig + 2Qi g Qi G ) 68 - 16Q Ig - Qig + 2Qi g Qi G 


(A32) 




80(1-Q G ii) 2 40(1 -Q G ii) 


P.3 
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56 - 13Qi g + 2Qi G + 38Qi g QiG 68 - 16Qi g - Qig + 2Qi g Qi G 


(A33) 




80(1-Q G ii) 2 40(1 -Q G ii) 


/3 




1 + 


35 - 4Qi g - Qig - 10Qi g Qi G 68 - 16Qi g - Qig + 2Qi g Qi G 


(A34) 




40 (1 - Qgii) 2 40 (1 - Qgh) 


S3 




1 + 


3 (16 + 5Qi G + 6Qi 2 G ) 31 - 4Q :G + 4Q 2 G 


(A35) 




80(1 -Qgii) 2 20(1 -Qgii) 
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1 + 


3 (27 - 17Qig + 6Qi 2 G ) 31 - 4Q :G + 4Q 2 G 


(A36) 




80(1 -Qgii) 2 20(1 -Qgii) 


«3 




1 + 


25 + 13Qig + 8Qi 2 G 31 - 4Q :G + 4Q 2 G 


(A37) 




40(1 -Qgii) 2 20(1 -Qgii) 


33 




1 + 


46 - 29Qig + 8Qi 2 G 31 - 4Q :G + 4Q 2 G 


(A38) 




40 (1 - Qgii) 2 20 (1 - Qgii) 


k. 




1 + 


34 - 17Qi g + HQ 2 g 37 - 13Qi g - 2Q 2 g 


(A39) 




40(1 -Qgii) 2 20 (1 - Qgii) 


h 




1 + 


37 - 5Qi g - 4Q 2 g 37 - 13Qi g - 2Q 2 g 


(A40) 




40(1 -Qgii) 2 20 (1 -Qgii) ' 



The coefficient 03 is achieved through the restrictions Xa^Xxp and XanX\%'i ^3 is achieved through XavXxp and XauXx^'i c 3 
through xijXffjl and x'xpX^ d:i through X^Xfit and x[$X^; e 3 through xLlxl$ and xilxl%; /3 through xi%xl$ and 
xi%xi°\ 93 through X^xf$ and XfiuXij^ h 3 through x-^x'Jp and X%xl°\ h through x^xi$ and Xplxj^; h through 
xt°Xpjl and x'Jpx'^ through Sp^xilXxly and & ltl ,xilx%\ and finally h through 5puX%x{% and 5~, v x%X%- 
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Finally, for those with three orientation dependent correlations, there exist 4 unique coefficients: 

i 74 - 9Qi g + 23Q? g - llQigQic - 20Q? g gi G 77 - 13Qi g - 2Q 2 g 

ai = 1 + 80(l-Q G n) 2 40(1-Q G „) (A41) 

1 74 - 20Qi g + 3Q : 2 g + llQigQic + 20g? g Qi G 77 - 13Qi g - 2Q 2 g . 

bi = 1 + 80(l-Q G n) 2 40(1-Q G „) (A42) 

75 - 12gig - 8Q 2 g + 2Qi G + 3Qi g Qi G + 10Q 2 g QiG 77 - 13Qi g - 2Q 2 g 

C4 " + 80(1-Qgii) 2 40(1 -Qgii) ( } 

i 77 - 9Qi g + 2Q 2 g - 2Qi G - 3Qi g Qi G - lOQfgQic 77 - 13Qi g - 2Q 2 g 

d4 = 1 + 80(1-Q G „) 2 40(1-Qgii) ■ (AM) 

The coefficient an is achieved through the restrictions XaLxlrfXy/jt an d Xa^xl^X^fii ^4 i s achieved through Xaix^x'Jj and 
xllX^x'tS; C4 through xi^X^X^ and X%X%X^pt\ and finally d 4 through X%X%xl° and x%X%xl%- 

For typical suppression ratios Qg// ~ 2/5, Qjg = Qig = 1/2, this produces the numerical values given in Table Al. 
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